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On partial derivatives 
of multivariate Bernstein polynomials* 
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Abstract 

It is shown that Bernstein polynomials for a multivariate function converge to 
this function along with partial derivatives provided that the latter derivatives 
exist and are continuous. This result may be useful in some issues of stochastic 
calculus. 


1 Introduction 

Widely known is the proof of the polynomial Weierstrass theorem based on Bernstein 
polynomials and on the Law of Large Numbers for Bernoulli trials proposed in [3]. 

This method is applicable for the approximation of multivariate functions, too. 
In the literature it was noted that for a univariate function, Bernstein polynomials 
approximate this function along with its derivatives assuming that they exist, see 
for example It turned out that this observation about an approximation 

along with derivatives holds true for multivariate functions as well; however, the 
authors were unable to find an exact reference and this was the reason for writing 
the present paper. In particular, this property is very useful in one of the proofs of 
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multi-dimensional Ito’s formula, which formula is the main tool in stochastic analysis 
(cf. [HI Theorem 10.4]). 

The latter proof starts with a verihcation of Ito’s formula for the products; then it 
is extended by induction from linear functions to arbitrary polynomials and finally to 
all functions with two continuous and bounded derivatives; this approach essentially 
follows the lines of one of the proofs of the Weierstrass theorem (cf. [T6l ch. XVI, 
§4]). It is essential that not only the function itself, but also its derivatives up to 
the second order were approximated on any compact. The latter requirement makes 
this version of the Weierstrass theorem a bit non-standard. At least, the majority of 
textbooks (cf. [161 [H]) where approximations by polynomials are discussed, usually 
focus on functions themselves and not on their derivatives. 

In the univariate case this nuance is not important because we can approximate 
the u-th order derivative and then integrate it n times. But in the multivariate 
case this trick does not help because integrals may depend on the paths. In [6] the 
property of approximation of a function along with its derivatives is simply claimed 
as widely known. The authors found it interesting to inspect whether celebrated 
Bernstein polynomials, which gave rise to a notable branch of approximation theory 
(cf. [H]) admit this property. The paper [TT] presents briefly the main result and the 
idea of approximation of partial derivatives in the case of d = 2. In this paper a full 
proof of this fact in the general case of d > 2 is provided. The case d = 1 is briefly 
presented in the next paragraph for the sake of completeness. 

This paper consists of hve sections. The hrst is Introduction. In the second 
some classical theorems and well-known generalizations in one-dimensional case are 
recalled. The third contains the main results about multivariate Bernstein polyno¬ 
mial derivatives convergence. The section 4 and 5 present the proofs of an auxiliary 
lemma and of the main results, correspondingly. 

2 The case d = 1 

For any function / on [0,1], approximating Bernstein polynomials are given by the 
formula 

:= 0 < I < !■ 

J=0 

Theorem 1 (Bernstein) // / G ^^([0,1]), then Bn{f',x) — )■ f{x), n ^ oo and this 
convergence is uniform [0,1]. 

See [3], and some generalisations in [12]. The following similar result holds true for 
the derivatives. 
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Theorem 2 /// G C^([0,1]), then Bn\f]x) n —)■ cxo and convergence 

is uniform on [0,1]. 

The proof of the Theorem [2] see, e.g., in 013; the result was established by 
I. N. Chlodovsky [1]. (The latter reference [1] gives only the title of his talk at the 
All-Unions Mathematical congress; the texts of most talks have not been published. 
According to m Chapter 1], since then the result of the Theorem [2] became known 
in the literature. However, the issue of correct references is a bit unclear.) There are 
various bounds of convergence rate under additional assumptions about smoothness 
or without them (cf. [13 (3 |TU] et ah), but they are not the goal of this paper. The 
following identity for derivative will be useful in the sequel (cf., eg., 0[Z])- 


n—1 


A -x) 

i=o ^ 


n-l-j 


( 1 ) 


where 

Later for a multi-dimensional case we will use a more detailed notation which 

emphasises that the increment corresponds to the variable x*; for a hxed z = 1/n we 
will use a short notation By induction the representation of the derivative of 

order k follows: 


I / ■ \ 

cUx'd -x)”---'. (3) 

where A^ is dehned also by induction as the operator A^ applied k times. For 
example for k = 2 we have, 

= A^{AJ{x)) = /(x + 2z) - 2/(x + z) + /(x). 

Note that the latter expression is one of the versions of a non-normalized hnite 
difference Laplacian for the one-dimensional case. Further, since n^A^^^f{x) =1 

7 -^ 1, n — )■ cxo, under any hxed /c, then by virtue of (I3l) and due 

[n — k)\ 

to the main calculus in the proof of the Theorem [1] based on the law of large numbers 
in the Bernoulli trials scheme with probability of success x (we recall that x G [0,1]), 
the proof of convergence for the derivatives follows immediately, of course, under 
the assumption / G C'^(R). It should be noted that all convergences are uniform 
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on [0,1]. We do not show here the details of this well-known calculation (cf. [HI 
Theorem 1.8.1] in the case /c = l) as for the main result of this paper in the multi¬ 
dimensional case they all will be given in the next sections. Yet we want to point 
out the development of the key Bernstein’s idea, which suggests to reformulate ([T]) 
and ([3]), respectively, as 

= ^nAi/nf 

and I 

where ^n{x) denotes a random variable with Binomial distribution Bin(n, x). These 
representations clearly explain why in the case d = 1 the left hand sides of ([T]) and 
(j3|) tend by the law of large numbers to their limits f'{x) and respectively, 

under the assumptions of the Theorem [2l 


3 The case d> 2 — Main Result 


Now let us consider d > 2, a; G M'^, x = {xi,, Xd). There are at least two ways to 
dehne Bernstein polynomials in the multi-dimensional case (actually, there are many 
ways and we will describe them later): either on a simplex := {[xi,X 2 ,... ,Xd) : 
0 < Xi,X 2 ,..., Xd, ||a;|| < 1) by the formula 


Bn{f;x) 


\ n 

0<Jl+i2H- 

ii 


h 

n 






(4) 


where j = (ji, ■■■ ,3d), and where a vector x norm, and the “modulus” of a multi¬ 
index j, and the polynomial coefficient (“n choose j”) are dehned as ||a;|| = 
\xi\ + \x 2 \^ -h \xd\, \j\ = ji + j 2 H- \-jd, and 

71 ^ 

r^j = _. 

ji\j 2 l---jd'-{n-\j\)\' 

or on a “d-dimensional square” (cube, etc.) by a formula similar to but different 
from (jll). 


Bn{f;x) 


Y / 

jl,j2,-,jd=0 



n 




X xf (1 - Xi)^-^^xi\l - X 2 )"-^^ ■ --xfil - Xd)”-'T (5) 
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Note that the degrees of the polynomials Bn and Bn are different. The formnlae 
(HI) and ([5]) allow probabilistic representation 

( 6 ) 

and 

Bnif] x) := Ef{n~^^n{x )), (7) 

where the distribntion of the random vector rin{x) of dimension d is the projection 
onto the hrst d coordinates of the polynomial (multinomial) distribution of dimension 
d+1 with parameters n (the nnmber of trials) and the vector of “snccess probabilities” 
{xi,... ,Xd, Xd+i), which satisfy conditions Xi > 0 {i = 1,... ,d+l) and J2i=i Xi = 1, 
while the random vector (,n{x) = {^h,---,^n) consists of independent components 
distribnted binomially Bin(n, Xi) each. 

Dne to the law of large nnmbers for Bernonlli trials for each coordinate we have, 

X, n -)■ oo, (8) 

n 

both in probability and almost snrely. For the seqnence {pn) the law of large nnmbers 
also holds trne: ^ 

—r}n{x) X, n ^ oo, (9) 

n 

in probability and almost snrely. The easiest way to show this is to use the law of 
large numbers for each coordinate, which follows directly from the one-dimensional 
version of this theorem for the binomial distribution. 

It apparently does not help too much to analyse partial derivatives where certain 
series still need to be treated. However, after the differentiation of these series we 
will get expressions, which admit representations via expectations of hnite differences 
for the function / with random vector arguments distributed either according to a 
polynomial law - or, more precisely, its projection - or as a direct product of binomial 
distributions, which will eventually lead to the desired result. 

Remark 1 Other variants are possible, such as a combination of a “cube” and of a 
simplex for different variables. For example, in case of d = di + d 2 and di > 2 one 
more version of multivariate Bernstein polynomials has a form. 
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Bn{f;x) := 




/ 


0<ii+i2H- 

ii j2vjdi>o 


JlvJdi 71 72 


- 

J di Jd^-\-3 

J d\ -\-d2 

\ n 

. . , , ? • • • 

n n 

n 


( di ji- 

Jdp 

•<‘| 

i-E-d 



V 7=1 / 


3 d^-\-\ 
^dl+l 

(1—Xi)"'“-^‘'i+^ • • • 

3di-\-d2 

d\-\-d2 


X a 


where x = (xi,... ,Xdi,Xdi+i,... non-negative, xi + • — hXd^ < 1, 

and 0 < Xdi+i < 1,. .., 0 < Xd^+dj < 1. 

It is quite likely that for such representations analogous results about convergence 
of polynomials and their derivatives may be established. 


So, a multivariate analogue of the Theorem [T] for the polynomials Bn and Bn may 
be formulated as follows. 


Theorem 3 /// G then 

Bn{f]x) /(x), X G 
Bnif-,x) ^ f{x), X G 

as n ^ oo. All convergences are uniform on and respectively. 

Note that any continuous function on or on can be extended to a continuous 
function on (of course, not uniquely). Convergence of both versions towards a 
continuous function / on the simplex or on the d-dimensional square/cube follows 
from the (multivariate) law of large numbers: in the second case it is applied to a 
sequence of independent and equally distributed random vectors with independent 
components, while in the hrst case - with dependent components. The proof can be 
easily found in may papers and textbooks (cf. e.g., [5]) and we omit it. 

To state the main result, let us introduce the following notations: : 

1. m — a natural number; 

2. k = {ki, k 2 ,... ,kd) — a multi-index; 

3. with a multi-index k k denotes the class of functions that possess a mixed 
partial derivative of order k = (fci, ^ 2 ,..., k^), which is continuous. 
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Recall that (C'™(M'^)) with m = 1, 2,... is a class of functions on with all 

well-defined mixed derivatives of order m, which are continuous (respectively, con¬ 
tinuous and bounded). The notation stands for the set of all mixed derivatives 
of the function / of order m. 

Theorem 4 1. If m > 0, f ^ then 

as n ^ oo, and all convergences are uniform on K'^ and respectively. 

2. If k = {ki, k 2 ,..., kd) is a multi-index and f G then 

as n ^ oo, and all convergences are uniform on K'^ and respectively. 

It is worth mentioning also the papers P and P devoted to asymptotic expan¬ 
sions and Taylor’s expansions for Bernstein polynomials of two and many variables. 
While being conceptually close, the results of the present paper do not follow di¬ 
rectly from these papers, and we do not aim to get neither asymptotic nor Taylor’s 
expansions here. 

4 Auxiliary result 

The following Lemma will be used in the proof of the main Theorem when a uniform 
convergence of hnite differences just under the condition of existence and continuity 
of the limit expression is needed. While elementary, this Lemma is required for the 
correctness of references and for a completeness of our presentation. 

Let {zi,... ,Zd) G M'’*, and (xi, ..., Xd) G M'’*. Let 

Ai/n,xJ{x) := f{xi, ..., Xi-i, Xi + -, Xi+i, ...,Xd) 

n 

f (^1) ■ ■ ■ 1 Xi—l, Xi, . . . , Xd), 

d f{xi, ...,Xd)-.= ■■■,Xd). 
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Lemma 1 For f G with k = {ki,..., ka) and x = {xi,... ,Xd) the following 

equality is valid 



For example, for ki = 0, integration over ^ 2 ) etc. is just not performed, and 

• • • > is treated as /(xi, ..., In particular, if all ki = 0 then the 

equality turns into the identity f{x) = f{x). 

The equality for two variables x = {xi,X 2 ) and /ci = /c 2 = 1 takes the form (for 
the sake of simplicity we omit the lower indices in 

rxi+zi l'X2+Z2 

Jxi Jx2 dxidx2 

The proof of the Lemma follows straightforward from the (one-dimensional) first 
theorem of the calculus (also known as Newton-Leibniz formula) by induction. 

5 Proof of the Theorem |4] 

0. It suffices to prove only the second part of the Theorem 0] because the first part 
follows immediately due to the identity 

= Pi C\R^). 

k: \k\<m 

1. Analogues of the formulae (1) and (2) in the multivariate case for Bn for the 



multi-index k = (/ci,..., /c^) could be written as 


^ 1^1 


Bn{f;x) = 


Q\k\ 


;Bn{f;x) 


dxk dxl\..dx^/ 

' n... {n — ki + 1) X n... {n — k 2 + 1) X ■ ■ ■ X n ... {n — kd + 1) 


X 


0<ji<n-ki, 


77,1^1 

ji h 3d 
n 


(jh qH . . . Qid 


X xf (1 - - X2T~'^^~^^ • • • - XdT~’^‘^~^\ 


( 10 ) 


where 


^ 1^1 


g\k\ 


dx^ ■ dx’l^ ... ex’"/ 




2=1 


and where ^{xg rneans A( 3 ..) = Ai/n^xi- 

The usage of the notation A^ with a multi-index k is correct because all operators 
A(a,.) commute. Indeed, let us denote (assume i < j for the sake of dehniteness) 

x^’^ = [x^, ... ,x^~^,x^ -f ... ,x^~^,x^ + 1/n,... ,x'^), 

f = [x^,.. . -I- 1/n, ... ,x'^), 


X 


X-' 


..., x^ -1- 1/n,..., x'^). 
For i ^ j we get elementary identities: 


A(^,)(A(^.)/(x)) = A(^.)(/(x^) - /(x)) 

= f{x"’^) - f{x") - {f{x^) - f{x)) 
= f{x"’^) - f{x") - f{x^) + f{x) 

and similarly, 

A(x7)(A(^,)/(x)) = A(^.)(/(x*) - /(x)) 

= f{x"’^) - f{x^) - {fix") - fix)) 
= fix"'^) - fix") - fix^) + fix). 
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Analogously to the representation (El), the equality flTOl) admits the following proba¬ 
bilistic meaning: 

91^1 ~ . ( n ... (n—ki+l) X n ... {n—k 2 +l) x ■ ■ ■ x n ... (n—kd+l) 

- 

X En"^lAV(n-'en-fc(x)); (11) 

the dehnition of the random vector ^n{x) = ..., see in §[3]). 

2. Let us prove the equality (1T0|) by induction. For A: = (0, • • • ,0) (a multi-index 
of order d) the desired equality (jS]) is equivalent to the dehnition of the polynomial 
this will serve as the basis of induction. It should be noted that for several 
variables (d in our case) induction can be carried in turn for the hrst variable up to 
ki, then for the second one up to ^ 2 , and so on. 

In other words, “double” induction (it can be also named multivariate) over each 
of the variables Xi, and then over the indices i = 1, ... ,d may be applied in our 
situation. For the induction step, it actually suffices to verify that the formula (1^ 
remains valid when any component of the multi-index k = {ki, ■ ■ ■ , kd) increases by 
one. For the sake of dehniteness let us check the step ki ki + 1 for i = 1. Recall 
that 


Jl J2 
n' 


Jd 


n 


nh nj2 r<jd 


BJf;x):= f 

jl,j2,-,jd=0 

X x{^ (1 - (1 - 


Because of a certain clumsiness of the formulae and for the sake of clarity we 
shall start with ki = 0 and then proceed to the general case. (The same approach 
will be used for the polynomials Bn). We get. 
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^^ \ n n n I 


^ Ji 32 3d 

jl,j2,---,jd=0 

3i 32 jd" 


E / 


0<ji<n, 

0<j2,...,jd<n 


(^ 3 in 32 . . . (^3d 

n n n 


X jixf ^(1—Xi)"' ( 1 — 0 : 2 )”' ^'^ ■ ■ ■ X^/{1—Xd)"' 

'jl 32 jd' 


njinn . .. ftjd 

, , . . . , I 

n n n 


- E / 

0<ji<n, 

0 <j 2 ,---,jd<n 

X x{^(n—Ji)(l—X2)”“'^^ • • ■ X^/{1—Xd)'^~^'^ 

Ef ( ^ I QJi+^Qj2 . . . ^id 

In’n’ ' n J ^ ” ” 

'll <r^'M '' 


0 <ji<n 

0 <j 2 ,---Jd<n 


X (jl + l)x{^(l—Xi)” X2)” • •X^‘^(l—Xd) 




E / 


(^3in32 . . . 

n n n ' 


0<ji<n, 

0<j2,---,jd<'n 

X x{^(n—ji)(l—xi)”“-^^“^x^^(l—X2)”“-^^ • • •x^‘^(l—Xrf)”“-^‘' 
Jl J2 3d' 




0 <ji<n, 

0<j2,--;jd<rt 


1 •>••••> 

n n n 


nn (^32 ... (^3d 
'^n-l'^n '^n 


xxi^(l—Xi)” ^ ■^^X2^(l—X 2 )” -^^ • • • x;^‘*(l—Xrf) 


>^-ld 


due to the identities 


(Ji + l)C?+‘ = - 


(ji + l)n! 


n\ 


(jl + l)!(n- (jl + 1))! ji!(n - 1 - jl)! 


= nCi\ 


n—1'! 


{n-33)Cij = - 


{n - ji)n\ 


n\ 


Ji!(n-ji)! ji!(n-l-ji)! 


= nCiU- 


Now, assuming that the formula ([3]) holds true for some k = {ki,..., kd), let 
us differentiate it once more in variable xi so as to get a similar formula for the 
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multi-index (fci -|- 1, /c 2 ,..., kd)- For the sake of brevity let us denote 

■ ■ ■ xj'(i - 

Note that since we are interested in the statement of the Theorem for n ^ oo, 
we can assume that all ki <n (although, for ki > n we have d^^Bn{x) = 0, and the 
right side of (|3]) also equals zero by dehnition of the number of combinations with 
negative n — kij. 

Derivative of a constant being equal to zero, we have. 
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/ 


d 

dxi 


\ 


0<ji<n-ki, 
\ i=l,...,d 


^ f ( ^ ^ I 


5: AV 

\ n n n 


d 

dxi 


\ 


. 0<ji<n-/ci 

y 2=1,...,d 




\n-ki-ji j 2 ,--;jd 

“fc 2 ,...,fcd 


(9 


n—fci 


axi \n n n J ^ 


jl=0 0<ji<n-ki 
i=2 . d 


\n-ki-ji j 2 ,--;jd 

“fc2,...,fcd 


n—k\ 




jl=l 0<ji<n-ki, 
i=2,...,d 

n—ki — 1 


\n-ki-ji^j2,--;jd 


k2,...,kd 


- E E 

il=0 0<ji<n-ki 
i=2,...,d 


kr Jlj2 Jd 


5 5 * • • 5 

n n n 


-fci 
n—fci —1 




Yj 1+1 J2 Jd 


E E ^‘/( = 


0<ji<n-ki, j'=0 

2=2,....d 


? ? • • • 5 

n n n 




J l'+l 
-fcl 

n—fci —1 / . . 

^ ^ \n n n 

0<ji<n-ki, ji=0 
i=2,...,d 


X C-i. i^‘("-*-'i-ii)(i-^i)“ ‘“S’Z'S, 


E 


«■ 


j2,---,Jd 

k 2 ,...,kd 


0<ji<n—ki^ 

1=2....,d 


\ \ n n n 

\ h=o ^ 


X 


''n—ki 




Jd 

n n n 


X CYfe, (n-/ci-ji)(l-xi) 


fi-fci-ii-i 


«■ 


J2,■■■,3d 


n—ki — 1 

0<ji<r2-/i:i, \ ji=0 

2=2,...,d 


E 


n-(fci+l)-ji 


X A 


"/[Zlll, _A‘/ 

\ n n n J \n n n 


E 


n—fci —1 


a 


.32,■■■,3d 

fc2,---,fcd 


^ (n-/ci) 


0<ji<n—ki 


ii=o 











as required. We have used the identities 

(n- fci)!(ji + 1) 


+ 1 ) - — 


(ji + l)!(n- fci - ji - 1)! 

^ (n-fci)(n-fci -1)! 

ji\{n — ki — ji — 1)\ ^ n-(fci+i)’ 

r<h ( u {n - ki - ji){n - ki)\ i ^ nh 

cu (" - - J.) = =(«-'=■) cL(..+i)' 

Hence, it follows by induction that the formula (E]) holds true. 

3. As for Bn, it remains to note that for / G due to the Lemma [H 

appropriately normalized hnite differences are uniformly close to the corresponding 
partial derivatives, that is, 

^ ^ I ix) X e K^, n ^ oo. 

^ ' dx'^ ^ ' 

So, the statement of the Theorem for the polynomials follows from flTT]) and 
from the law of large numbers ([8]) . 

4. The analogue of the formula (jS]) for the polynomials Bn has a form. 


Q\k\ Q\k\ 

Bi^\f-x) = —Bn(f]x) = -r-r- 

Q^k dx’l^dx^^\..dxf 




E A^l A ^2 \^d 

■■■<n—\k\ ^ 11/ 


0 <ji+j 2 +---<n-\k\ 


P ; Jl J2 Jd, \ rJi rJd(-\ 

^ I LT’ W’ • • • ’ W ^n-\k\ ^1^2 "'^d 


3d 


JI n nj 
For the seequel, let us recall the following notation. 


\x 


n-|fc|-|j| 


( 12 ) 


j\ki \k2 /\kd _ Afc 

^ ■ 

The formula flT^ also has a simple probabilistic meaning. Similarly to the rep¬ 
resentation m we write. 


B<‘>(/;n = E- 


n! 


{n-\k\) 

for the dehnition of a random vector rin{x) see §[31 


A^/(n ^r]n-k{x))-, 


(13) 
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5. For the proof by induction of the relations flT^ and flT^ . let us note that the 
basis of induction {k = 0) coincides with the dehnition of the polynomial Bn{f]x)] 
as for the inductive step, for the sake of clarity of the following bulky computations 
let us hrst differentiate once the function Bn{f]x). we have, 


d^^Bn{f-,x) := f 


0 <jl+j 2 -\ - \-jd<n, 

ii j2vid>o 


n J2 3d 

1 ? • • • 5 

n n n 


X Cldx^X^xf ■ ■ ■ X^^{l-Xi-X2 - XdY 


o<ii+i2H— \-jd<n, 
ji>o, i2v,id>o 


^ ' n n n 


X • • ■ X’^{l-Xi-X2 - Xd) 




E / 


0 <il+i 2 H- \-jd<n, 

jl,j2:---,jd>0 


3j_ 3JI j_d 

1 5 • • • ? 

n n n 


F't? • • • 'r^3(^ —Ti —Tr.— . . . — Xd}^ ^ 


X Cl[n—\j\)x\^xf ■ ■ ■ x^/{l—xi—X2 

E 


0 <il+t 2 H- \-jd<n-l 

jl,j2,---Jd>0 


./jl + 1 j2 jd 

J \ ? 5 • • • 5 

\ n n n 


—Ti —'To— . . . — Xd)^ 


X X2 ■ ■ ■ x^^{l—xi—X2 


0<il+i2H- \-jd<n, 

jl,j2,---,jd>0 


E / 

^ \ n n n 


• • • 'rB(^ —'Ti —'To— . . . — Xd}^ ^ 


X (n—\j\'jx\^X 2 ^ ■ ■ ■ X^/{1—Xi—X2 


0 <ji+j 2 -\ - \-jd<n, 

ii J2,...Jd>0 


A r Ml -^2 jd 
\n n 'toy 

X Cl^3^i’^‘^x{^xi^ ■ ■ ■ xjf{l—Xi—X2 - 
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due to the identities 


+ 1 ) 


^KJi +1) ^ ^ (n-1)! ^ ^n,n,...,u 


Cl 


1 J2,.--Jd 


{n-ji - jd) 


n\{n-ji - jd) ^ 


6. The “full” induction step: under the assumption that the formula flT^ holds 
true for some multi-index k = {ki,... ,kd), let us differentiate this expression again, 
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say, with respect to xi. We get, 


d„Bf\f-,x) = 4. 




0 <ji+j 2 -\ - l-id<n-|fc| 

^ j'2, ■ ■ ■ ,jd ™i2 (1 


— — — 
n — |/c|)! \ n ’ n ’ ’ n 


X 


n-|fc|-|j| 


E 


n\ 


0 <ji+j 2 -\ —l-id<n, 
jl>0,j2,--;jd>0 




n 


J2 


n n 


Jd 

n 


X C^n-\k\j^^l ^^2 • • • X^/{l-Xl-X2 - Xd) 


n-\k\-\j\ 


E 


n\ 


0 <il+i 2 H- \-jd<n, 

Jl J 2 ,...Jd >0 


{n-\k\)\ 




n n 


Jd 

n 


X C^n-\k\{'^~\^\~\j\)^i^^2' ■ ■Xd^il-Xi-X2 - Xd) 

E 7-4r^AV 


i-lil-1 


0<il+i2H-l-id<n- 

jl,j2,---jd>0 


, (n-|fc|)! E' 


fc ^ Jl + 1 J2 jd 

? ? • • • ? 

n n n 




E 


ri! 


ks-l ji J2 3d 


0 <ii+i 2 H-hid<n, 

jl,j2,---jd>0 


{n-\k\)\ 


AV 


) ) 
n n 


n 


X xf • • •a:^''(l-xi-X2- Xd) 

.5:. 

0 <J 1 +J 2 H- \-Jd<n, VII/ ^ 

ii,i2,...id>o 

E 


i-|fc|-i-|j| 


n 


0<ii+j2H- ^jd<n, 

jl,j2,---jd>0 


X X 2 x^ {1—X1—X2 

by virtue of the identities 


( I 7 I H '1 I ^( 3 : 1 )'^ / I ) 1 • • • 1 

[n — \k\ — 1)1 \n n n 


—T,—'r,7— •••— Xd)^ ^ 
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'^n-\k\ 


(Jl + 1) = 


{n- \k\)\{ji + 1) 


= [n 


(ji + l)!j2!---(n - \k\ - \ j\ -l)! 

{n — \k\ — 1)! 






ji'jj!--' (n - |fc| - 1 - |j|)! 

_ (n - |fc|)!(n- |fc| - |j|) 

,jj{n-\k\-\j\y. 


Hence, by induction the formula (IT^ holds true along with (IT^ . 

7. Again having in mind how Bernstein’s method for the Theorem [T] was applied 
in the proof of the Theorem 0] for the polynomials Bn, and similarly to what was 
done earlier for Bn, let us recall that for / G the normalized hnite differences 

are uniformly close to the corresponding partial derivatives, namely, 

^ ^(U. for xesy „ ^ oo 

(cf. the Lemma [T]). Hence statement of the Theorem for the polynomials Bn follows 
from the law of large numbers (jD]) and from the representation flT^ . 

The Theorem |4] is proved 
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